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1.  Introduction 


This  report  provides  all  of  the  necessary  details  of  the  analysis  of  dipole  antennas  printed  on 
grounded  magnetodielectric  substrates.  Since  the  sources  considered  are  essentially  impulses, 
suitable  Green’s  functions  could  be  defined  from  the  formulas  contained  to  facilitate  analysis  of 
any  radiating  sources  by  superposition.  The  interest  in  magnetodielectric  substrates  stems  from 
the  expectation  that  new  materials  possessing  specially  tailored  magnetic  properties  may 
revolutionize  the  engineering  of  thin  planar  antennas.  Since  these  materials  often  require 
complicated  constitutive  equations  to  describe  their  properties  rigorously,  the  solution  of  intricate 
boundary  value  problems  has  become  important  to  the  understanding  of  the  fundamental 
properties  of  associated  radiating  structures. 

The  boundary  value  problem  considered  in  this  report  is  simplified  substantially  by  assuming 
that  the  magnetodielectric  may  be  modeled  as  a  simple  material  with  scalar  constitutive 
parameters.  This  approximation  is  likely  not  entirely  satisfactory;  however,  it  does  reduce  the 
problem  under  consideration  to  a  well-known  boundary  value  problem  in  applied 
electromagnetics.  Despite  the  extensive  amount  of  existing  literature  treating  this  problem,  no 
completely  self-contained  treatment  of  the  mathematical  boundary  value  problem  itself  was 
found.  For  this  reason,  this  problem  is  given  a  full  treatment  to  demonstrate  the  method  of  attack, 
which  may,  in  future  research,  be  extended  to  treat  the  case  of  magnetic  substrates  described  by  a 
more  realistic  set  of  constitutive  equations. 

Following  Arnold  Sommerfeld,  we  pose  the  boundary  value  problem  in  terms  of  magnetic  vector 
potential  and  proceed  to  solve  for  the  unknown  coefficients  by  applying  the  boundary  conditions 
(the  complete  derivation  of  which  may  be  found  in  Appendix  A).  In  the  second  section,  the 
electromagnetic  field  is  calculated  from  the  magnetic  vector  potential.  Here,  we  expend  some 
additional  effort  to  reduce  the  form  of  solution  to  an  especially  compact  and  convenient  form. 
These  solutions  may  be  compared  for  accuracy  to  the  solutions  reported  in  the  references. 
Sections  3-6  are  devoted  to  further  analysis  and  manipulation  of  the  solutions  to  obtain  results 
useful  in  practice  such  as  surface  wave  excitation,  far-field  radiation,  directivity,  radiation 
resistance,  and  efficiency.  Finally,  it  should  be  noted  that  a  novel  approach  is  used  in  the 
treatment  of  the  complex-valued  square  root  functions  in  the  theory  of  complex  integration  used 
to  derive  the  surface  waves  excited  on  the  magnetodielectric  substrate.  Appendix  B  shows  how  a 
conventional  definition  for  the  square  roots  can  be  given,  which  can  be  shown  to  provide  all  of 
the  necessary  (and  sufficient)  properties  to  enable  a  contour  to  be  closed  in  the  lower  half 
complex  plane. 
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2.  Solution  of  the  Boundary  Value  Problem 


The  boundary  value  problem  is  illustrated  in  Fig.  1 . 

A^o  ^  ^0 

z  =  0 - - - > - 

li,e 

^  //////////////// 


Fig.  1  The  boundary  value  problem 


z 


Here  we  have  an  x-directed  Hertzian  dipole  situated  at  an  air-dielectric  interface.  The  dielectric 
layer  has  a  thickness  h  and  is  backed  by  a  perfect  electric  conductor  (PEC)  plane  located  at 
z  =  -h.  The  current  density  associated  with  the  Hertzian  dipole  is 

/  =  mod(r)  (1) 

The  electromagnetic  field  generated  by  this  current  can  be  represented  by  2  components  of 
vector  magnetic  potential:^ 

A  =  -\-  zA^  (2) 

which  satisfies  the  homogeneous  Helmholtz  equation  in  both  regions,  but  not  at  the  z=0 
interface: 

(V^  +  k^)A  =  0  (3) 

A  general  solution  to  Eq.  3  may  be  obtained  for  the  Cartesian  components  of  A: 

00 

Ay^  —  jj  Aye~-'^  e~-'^yydkxdky 

—  00 

/Cz  =  Jk^  -kx  -ky 

The  solution  to  the  problem  consists  of  finding  the  coefficients  (a,  b)  that  satisfy  all  the  boundary 
conditions.  The  geometry  of  the  current  problem  suggests  that  we  set 

A^^  =  ,  z  >  0 

,  z  >  0  (5) 

Ax^  =  -I-  ,  —  h  <  z  <  0 
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A  =  fe 


h  <  z  <  0. 


In  order  to  determine  all  of  the  unknown  coefficients  in  Eq.  5,  we  require  6  boundary  conditions. 
The  derivation  of  these  conditions  is  somewhat  involved.  They  are  listed  below  (see  Appendix  A 
for  details). 


dA, 


dz 


^Zi 

z=0' 

A 


dA 


Zi 


dz 


^Xo\  ^  =  0 

^2  1  z=-h 

(6) 

o 

II 

1 

II 

N 

\ 

N 

(7) 

1  ^  ~  1 

1  z=0'*'  ~  ^  ■^Z2  1  z=0~ 

(8) 

1  saT  , 

. - z-0-  =  -ko£Io 

g-r  dz 

(9) 

M  ' 

N 

II 

O 

+ 

II 

( 

N 

II 

O 

1 

(10) 

dAz,  1  UyEy  —  1  -  . 

^  —ik  ^ 

dz  =  ^  ^  gy£y  = 

(11) 

Using  the  expression  for  A^  in  Eq.  6,  we  find 


Similarly,  from  Eq.  7, 


+  de  =  Q 
d  =  — 


g  = 


Using  Eqs.  12  and  13,  we  can  return  to  Eq.  5  and  update  our  expressions  as  follows: 

Ar.  =  ,  z  >  0 

z  >  0 

,  —  h  <  z  <  0 

=  f{e~j’^^2^  +  e^j’^^2’^ej’^^2^)  ,  -  h  <  z  <  0. 

Next,  we  plug  the  updated  expressions  Eq.  14  into  the  boundary  condition  Eq.  9: 

1 

-jk^^a - -jk^^e^J’^^2’^)  =  -go^lo- 

kr 


Az^  =  be 
=  c{e-^^^2^  - 


Similarly,  from  Eq.  10, 

Substituting  Eq.  16  into  Eq.  15, 


a 


c(l  - 


(12) 


(13) 


(14) 


(15) 


(16) 
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—jkz^c{l  —  H - ikz^c{l  +  =  —^LqIIq 

kr 

-;A:2^ce^^2z^(^e“^'^z2^  -  +  J_y/f^^ce;fcz2^(^e-;'cz2^  +  e^'^zz^j  = 

kr 

2 

—2kz^ce^^^'2-^  s\n{kz^li)  H - cos^^kz^h)  =  —fJ-o^Io 

kr 

2jcsm{kz^h)  e^'^^2f^[firjkz^  +  kz^  cot{kz^h)]  =  -fiofirflo 
We  denote  the  bracketed  expression  in  Eq.  17  as  follows:^ 


Dte  =  krJkz^  +  Kz  cot{kz^h) 


Then, 


c  = 


-liQlirilo 


2jDYE  sin(/c22^)  e^^^z^ 


Substitution  of  Eq.  19  into  Eq.  16  gives 


a  = 


-kokr^lo 


.(^1  _  e'^jkzzh^ 


a  = 


2] Dte  s\n{kz^h) 

-kokr^lo  ^jkz^h^^-jkzzh  _  ^jkzzh^ 


2jDTE  s\n{kz^h^  gl^zz^ 


a  = 


fiokr^lo 


D 


TE 


Next,  we  apply  boundary  condition  Eq.  8,  obtaining 


(17) 

(18) 

(19) 


(20) 


h  =  — /(l  +  e^^'^^2^) 

kr 


and,  from  boundary  condition  Eq.  1 1 

1 


-— /(-y^zz  ^jk, 


—  1 

£j^ 


■a 


(21) 


Substitution  of  Eqs.  20  and  21  into  Eq.  22  gives 

-yt..  2/(1  + 

f^y  Cy  j^y  Cy  U 'f  J7 

1  .  .  1  .  . 

“7^zi  — fe^^^2^{e  ^^^2^  +  +  jkz^ - fe^^^2^{e  ^^^2^  — 

j^y  jA,y£y 

[iyEy  1 /Tq/Tj-Z^/q 

=  jkx - j: - 

kr^r  kljE 

—2jkz  — /e-'^z^  cos(/c2  /i)  +  2/C2  - /e-'^z^  sin(/c2  /i)  =  - k^k21 — 

^  fir  ^  ^  ^^lirSr  ^  kr^r  D 

-2f  cos{kz^h)  e>^^2^  ^  h)]  -  ik  ~  ^  kokr^k 

\^£yJK.Zi  ^Z2  J 


(22) 


TE 

(23) 


IJ,y  £y 


pLyEy  D 


TE 
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We  denote  the  bracketed  expression  in  Eq.  23  as  follows:^ 

Dtm  ==  -  /fz2 


Then, 


/ 


2Dte  DTMCOs{k^^h) 


(24) 


(25) 


Substitution  of  Eq.  25  into  Eq.  21  gives 

1  flQfl^£lQjkj((^fl^£^  1) 


b  = 


b  = 


kr  2Dte  OrMCOs^kz^h) 


(1  + 


2DrE  DjMCOs{k^^h)  e^'^^2’ 


b  = 


~ko^^oj^x^kr^r  1) 


DjeD 


TM 


(26) 


This  completes  the  solution  for  the  unknown  coefficients.  We  can  now  return  to  Eq.  14  and  write 

-kokr^h 


^  — _ 

2jDTEsm{k^^h)ej’^^2>^ 


_  -kokr^lp 

2jDTE 


.^jkz^h^^-jkz^iz+h)  _  ^jkz^iz+h)^ 


^X2 


/ro/Tr^/o  sin  [/Cz^(z  +  /i)] 
Dte  sm{k^^h) 


(27) 


Similarly, 


kokr^Ipjkxikr^r  1)  r  -jk.^z  ,  2jk22h  jk22z\ 
2DrEDrMCOs{k,2h)e^’^^^’^^  ^ 

~  _  -kpkr^ipjkxikr^r  "  1)  COs[A:z^(z  +  /l)] 

Dte  k^TM  ^os  {kz2 


(28) 
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Assembling  all  of  our  results,  the  complete  solution  to  the  boundary  value  problem  is 


Dte 

>  0 

DteDtm 

z  >  0 

^^2  -  n  ■  fir  ,>,sin[/c2  (z  +  /i)]  , 

Dte  sm{k^^h) 

—  h  <  z  <  0 

T—  ~kokr^toj^x(^kr^r  ~  1)  t;  r  i 

A,  =  - 7 - r - COS  /C,  (Z  +  h)\ 

Dte  Dtm^^^{}^z2^ 

,  —  h  <  z  <  0 

Dte  =  krJkz^  +  cot{k^^h) 
Dtm  =  ^rjkz^  -  tan{k^^h) 

^Zl  ~  ->1^°  ~  kx  ~  ky 

kz2  ^  kr^rko  kx  ky 

(29) 


Note  that  a  2-dimensional  inverse  Fourier  transform  is  required  to  recover  the  spatial  dependence 
of  the  fields.  It  is  also  worth  mentioning  here  that  Djg  and  Djf^  represent  the  dispersion 
equations  for  supported  surface  wave  modes  that  are  transverse  electric  (TE)  and  transverse 
magnetic  (TM),  respectively,  to  the  z  axis.  As  such,  our  next  step  is  to  derive  the  z-components 
of  the  electric  and  magnetic  fields.  Most  quantities  of  interest  can  be  calculated  directly  from 
these  fields. 


3.  Z-directed  Electric  and  Magnetic  Fields 


The  z-components  of  the  electric  and  magnetic  fields  are  related  to  the  magnetic  vector  potential 
as  follows: 


Hz  = 


jky 


i4'v 


](jl>IJ.£ 


-JK 


dz 


+ 


(30) 


The  equation  for  the  z-directed  magnetic  field  is  rather  simple.  Using  our  results  from  Eq.  29,  we 
obtain 
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D 


-jk^^z  z>  0 


TE 


H. 


ilojky 


Z2 


DrEsm{k^^h) 


sin 


[k^^^z  +  h)]  ,  —h<z<0 


Next,  focusing  on  the  z-directed  electric  field  for  z  >  0,  we  first  write 


(31) 


-JK 


dA, 


dz  Dee 

Substitution  of  Eq.  32  and  the  expression  for  i42  from  Eq.  29  into  Eq.  30  yields 


(32) 


KqUte^TM 


(33) 


Although  it  is  difficult  to  tell,  the  form  of  Eq.  33  can  actually  be  considerably  simplified.  In  a 
(non-obvious)  step, 

(.kr^r  ^y  )  kr^rij^O  )  (^kr^r^O  ^Z2  ) 

+  ky^)  =  kz2^  -  kr^rK^^  (34) 

We  substitute  Eq.  34  into  Eq.  33: 

^loVojk 


^  ^oDeeD 


e  ^  zi  |y(^/c22  kr^rkz^  )  "f  krkz^^EM^ 

TM 

^loVojkx  _n 


■e  \jkz2  krk-z-^ij^EM  ^rJ^Zi)] 


^  koDrE^^TM 

E..  =  +  krkzX^rjkz^  -  kz^  tan{k^^h)  -  SrjkzJ] 


"Zl 


koDEED 


TM 


E..  -  firk^^k^^tan{k^^h)] 


"Zi 


koDEE^TM 
^^oVokxkz2 


^  koDEED 


"Zi 


tan{k^^h)  e  cot{k^^h)  +  /r^y/CzJ 

TM 

^^oVokxkz 
koDEE^ 


'■Z2 


TM 


tan{k^^h)e  ^’"^^^[Dee] 


Ez,= 


floTjokxk 


Z2 


koD 


TM 


e  ,  z  >  0 


(35) 
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Lastly,  we  must  find  the  z-directed  electric  field  for  —h  <  z  <  0.  As  before,  we  first  write 


-iK 


dA 


X2 


-fiolir^Iojkxk 


Z2 


cos[k2^(z  +  /i)] 


dz  DTEsm{k^^h) 

Substitution  of  Eq.  36  and  the  expression  for  A^^from  Eq.  29  into  Eq.  30  yields 


(36) 


_  ^iQ^r^Ipjkx  r 

upon  cos[A:22U  + /iJJ 


XflrSr  -  l)ikx^  +  ky^) 


+  ■ 


^•22 


DjM^os{k2:^h)  s\n{k2:^h) 


(37) 


As  it  was  for  ,  we  can  also  simplify  Eq.  37  considerably: 


_  -a^kx  cos[/c^^(z  +  /i)] 

^  (xiEQEyDjpDjj^  cos(^k2^h^ 


[(Er^r  ky  )  +  k)TMkz2 


_  _  -ilpkx  cos[/c^^(z  +  /i)] 
^  ojeqE^DeeDyj^  cos(^k2^h^ 


[(^Z2  kr^rkzi  ) 


E..  = 


+  {^rjkz^  -  kz^  tan{k^^h))k^^  cot{k^^h)] 
-£Iokx  coslkz^iz  +  h)] 


(x)EqEj.DjeDjj^  cos{k2^hL) 

^  _  -ilpkx  cos[/c^^(z  +  E)] 
a)£o  COS^^k^^h^ 


\kz2  kr^rkz^  ”f  ^rj kz^k^^  cot^^/c^^h.^  k^^  ] 


[^r^r^Zi  ”f  ^rj kz^k^^  C0t(^/C22^)] 


-Erilokxk^^  cos[k^^{z  +  h)] 


E  = 

ojEqE^DeeDej^  cos{k2;^h) 

~  _  -flpjkxkz^  cos[/c^^(z  +  E)] 
ojEqDjeDjj^^  005(^22/1) 


[fipkz^  +  jkz2  cot{k^^h)] 


p-rjkz^  +  /fz2  cot{k^^h)] 


~  _  -^Ipjkxkz^  cos[/c2^(z  +  E)] 
^  (dEqDyeDj]^  cos(^k2^h^ 


Wte] 


Ez2- 


^IpVpjkxkz^  cos[/c2^(z  +  h)] 
kpD-TM  cosik^^h) 


-h  <  z  <  0 


(38) 


8 


Assembling  these  results,  we  get  the  following: 


_  HoUrjky 

Zi  n 

Ui 


'TE 


ilojky 


DTEsm{k^^h) 


sin  [/c^^Cz  +  h)] 


^loVoiKK^  cos[/Cz^(z  +  h)] 
^oDym  cosik^^h) 


z  >  0 

—  h  <  z  <  0 

,  z  >  0 

—  h  <  z  <  0 


(39) 


As  mentioned  previously,  in  order  to  determine  the  aetual  spatial  dependence  of  the  electric  and 
magnetic  fields,  we  must  perform  the  inverse  2-dimensional  Fourier  transform  of  the  expressions 
in  Eq.  39.  Unfortunately,  it  appears  to  be  impossible  to  integrate  Eq.  39  in  closed  form.  However, 
far  from  the  Hertzian  dipole,  we  can  employ  the  stationary  phase  approximation  to  the  integrands 
to  obtain  the  far- fields.  Additionally,  it  is  also  possible  to  extract  TE  and  TM  surface  wave  fields 
by  interpreting  the  integral  as  a  contour  integral  in  the  complex  k^  and  ky  domains. 


4.  Far-Field  Radiation 


Erom  Eq.  39,  we  could,  for  instance,  write  the  solution  for  the  z-directed  electric  field  above  the 
substrate  as 


00 

~  (2nY  II  e  j'^y^dky-dky 

—  CO 


(2^  II  tan{k,^h)  e  e  i^yydk^dky 


(40) 


The  integral  Eq.  40  is  of  a  general  form  for  which  the  stationary  phase  approximation  is  well 
known:^ 


where 


cos  6 


"Zi 


^qD 


-tan 


TM. 


{kz2,h) 


>-Ao?’ 

2nr 


,  (far  field) 


(41) 


9 


(42) 


Similarly, 


Zi 


II  j’^yydk^dky 


II  e-^-^z,Zg-;7c.x  e-Jhydk.dky 


(43) 


and 


H. 


ilolifky  COS  0  e 


Zl 


D 


TE, 


2nr 


Additionally,  under  the  stationary  phase  approximation. 


Er  =  0 


(far  field) 


(44) 


Eo  =  -■ 


sin  9 

u  _^0^z 
^d)  .  ^ 

^  sm  0 


(45) 


Therefore,  a  good  approximation  to  the  radiated  field  in  the  far- field  in  the  upper  hemisphere  is 


Er^O 

^JoVojkxskz2,cote  . 

Eg  = - — - tan  [k^^^hj 


k^jD 


TM. 


2nr 


Ecf,  = 


tlorjQiiyky  cot  0  e 


D 


TE. 


2nr 


(46) 
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5.  Surface  Wave  Fields 


We  now  turn  toward  the  derivation  of  the  surface  wave  fields  generated  by  the  Hertzian  dipole. 
In  order  to  do  this,  we  must  interpret  our  integrals  as  contour  integrals  in  the  complex  and  ky 

domains.  For  convenience,  the  integrals  are  repeated  below: 


(271)2 


aoM'ky  -jk^x  e-j'^yydk^dky 


1  rr  ^ipjky 

(27r)2  JJ  Dte  sin(A:22/i) 


sin[k22(z  +  h)]  e  e  ^^y^dk^dky 


(27r)2 


ff  ^  tan(/Cz^/i)  e  e  ^^y^dk^dky 

JJ  koDji^ 


(271) 


00 

^// 


^^oVojkrkz.  coslkzAz  +  h)] 

,  ^ e->^yydk^dk^ 

koDjM  cos{k,^h)  y 


We  introduce  the  change  of  variables, 

kx  —  kf-  cos  a 
ky  —  kf  sin  a 

(48) 

X  —  p  cos  0 
y  =  p  sin  0 

Noting  that 

kxX  +  kyy  =  /ctp(cos  a  cos  0  +  sin  a  sin  0)  =  kf-p  cos(a  —  0)  (49) 
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the  integrals  in  Eq.  47  become 

2n  00 


H. 


Zi 


=  C^// 


ilQfirjkt  sin  a 


H. 


Z2 


=  (5^// 


Dte 

0  0 

2n  00 

ilojkf-  sin  a 


^-jk,^z^-jktpcos{a-4>)  ktdktda 


0  0 
271  00 


"Zi 


=  -c^// 


sin{kz^h) 
iloTjokt  cosak^^ 


sm[k^^(z  +  /i)]  e-7'^tPcos(a-0)  k^dk^da 


(27r)2 


0  0 
2n  00 

// 

0  0 


koD 


TM 


tan^kz^h^  e  i^tpcos{a  <P)  i^^^dkf-da 


^^oVoj^t  cos  a  kz^  cos[kz^(z  +  /i)] 
hDrM  cos{kz^h) 


e-)'^tPcos{a-ct>)  j^^dktda 


(50) 


With  the  aid  of  the  following  2  mathematical  identities, 

2n 

— ZnjJ^ikf-p)  cos  <p  =  j  cos  a  da 

0 

271 

— 2717/1  (/Cfp)  sin  0  =  f  sin  a  da 


0 

the  integral  expressions  in  Eq.  50  reduce  to 


(51) 


(52) 
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With  the  integrals  put  into  the  above  form,  it  is  possible  to  close  an  expanding  semicircular 
contour  around  the  poles  (zeros  of  Dje  or  Dj^)  in  the  lower  complex  kf-  plane.  Complications 
arise  in  the  proper  definition  of  : 

which  cannot  be  single-valued  and  continuous  throughout  the  entire  complex  kf-  plane.  The  lines 
across  which  this  function  is  discontinuous,  the  branch  cuts,  depend  on  the  choice  of  the  positive 
or  negative  root  in  various  regions  of  the  complex  kf-  plane.  The  algebraic  form  of  Sommerfeld’s 
choice  for  k^^  is  discussed  in  Appendix  B  and  is  given  by 

=  -J|/co^  -kt^\exp  ^argo(/co^  -  (56) 

where  argo(A:o^  “  is  the  value  of  arg(/co^  —  /ct^)  li^s  between  0  and  2n.  The  branch 
cuts  and  behavior  of  k^^  for  this  choice  are  depicted  in  the  complex  kf-  plane  in  Fig.  2. 
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Re{kzJ  <  0 
lm(kzj  <  0 

Re(kzJ  >  0 
\m{kzj  <  0 

/Cq 

— /Cq 

Re(fczJ  >  0 
lm{kzj  <  0 

Re(kzJ  <  0 
\m{kzj  <  0 

Branch  Cut  — 

Fig.  2  Branch  cuts  and  behavior  of  for  this  choice  in 
the  complex  plane 


On  the  real  axis,  just  above  the  branch  cuts,  this  definition  reduces  to  the  usual  requirement  that 
waves  be  attenuated  exponentially: 


kz^  —  kf  ,  kf-  <  ko 

h 

Throughout  the  rest  of  the  complex  plane,  Eq.  56  preserves  the  property  that  Im(/c2^)  <  0.  The 
closed  contour  that  we  use  is  shown  in  Fig.  3. 


Fig.  3  The  closed  contour  that  we  use 
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By  the  residue  theorem,  we  have 

+  ^2  +  4  +  ^4  +  ^5  +  ^6  +  ^7  +  4  =  —2nj  ^  res(/Cp)  (57) 

Since  the  semicircular  arc  is  expanding  into  a  proper  region,  Jordan’s  Lemma  says  that  this 
integral  will  go  to  0  as  the  radius  of  this  arc  goes  to  infinity: 

^2+4^0  (58) 

Additionally,  the  integration  around  the  branch  point,  1^,  vanishes.  Therefore, 

k  =  -ih  +  h)  -  ik  +  k)  -  27r;  ^  res(/Cp)  (59) 

where  +  h)  ^nd  (1 4  +  4)  do  not  vanish,  since  is  discontinuous  across  the  branch  cuts. 
An  asymptotic  analysis  of  the  integration  around  the  branch  cut  by  the  method  of  steepest 
descent  gives  a  result  for  the  radiated  field  that,  under  certain  conditions,  is  in  agreement  with 
our  earlier  results  obtained  by  the  method  of  stationary  phase.  Thus,  Eq.  59  may  be  interpreted 
physically  as  a  decomposition  of  the  field  into  radiated  and  surface  wave  components. 


In  general,  the  locations  of  the  zeros  of  either  or  Dtm  must  be  found  by  a  root- finding 
algorithm  on  a  computer.  However,  analysis  of  their  expressions  indicates  that  our  contour  will 
capture  these  points  in  quadrant  IV  of  the  complex  kf-  plane.  Furthermore,  the  real  part  of  these 
zeros  will  be  located  somewhere  in  the  interval  (/cq  ,  V/Ur^r^o)-  lu  the  limit  of  a  lossless 
substrate,  /Tj.  and  Sj-  are  real  numbers,  and  the  zeros  migrate  up  to  the  real  axis  of  the  kf-  plane.  In 
this  limit,  it  should  actually  be  understood  that  the  zero  is  “slightly”  below  the  real  axis,  so  that 
the  formula  Eq.  59  still  applies.  Otherwise,  the  formula  Eq.  59  would  be  a  discontinuous  function 
of  the  substrate  properties  in  passing  from  a  small  amount  of  loss  to  no  loss. 

Using  the  expression  for  the  form  of  a  TE^  surface  wave  above  the  substrate  is  found  by 
calculating  the  residue  at  any  zero  of  D-pg : 


H 


sw  _  yA  _  J  J  vl  Tip 


Zl 


■e-'4  lim  kf-  —  ky 


Iht^-jk^yz^-jktpu 


/Cf  ~>/Cp 


D 


TE 


H. 


sw  iloPrSinip  jn  2kp 


Zl 


.^■>4 


np 


e  i^pPkr,  lim 


k-t^kp  Dpg{kk)  k)TE{kp) 


cw  £Ir,UrSin(b  /TT  2/Cr,  1 

sw  ^_0Mr_ - _^^-jk,yZ^-jkpPi^  - ^ 

^Drg{k^) 


(60) 
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where  Djg  (/Cp)  is  the  derivative  of  with  respect  to  evaluated  at  the  pole  location  kp.  k^^ 
is  also  evaluated  at  the  pole  location.  By  the  same  process,  we  can  write  the  expressions  for  the 
z-components  of  the  surface  wave  fields: 


TE^  Surface  Wave 


„  sw  iloErSincp 


np 


5w  ^/osin0 

Z2  9 


2/CpSin[/c,^(z  +  /i)]^_.^^^^ 


Dte  (^p) 
1 


np  sin^k^^h) 


Dte  (^p) 


TM^  Surface  Wave 


"Zl 


2/co 


Trp 


Dtm  (^p) 


sw 


"Z2 


noVoCos(p  jE  \2kp  cos[k^(z  +  h)]  _  1 

= - :r; - ^-'4  j - — e  ^  P^k^^kp - - 


2/co 


Tip  cos{kz^h) 


Dtm 


(61) 


6.  Radiated  Power,  Directivity,  and  Efficiency 


The  Poynting  vector  describing  the  power  flow  associated  with  the  far-field  radiation  is  given  by 
the  following  formula: 


1  1 

Savg  =  —ReiE-E*)  =  f-{\Eg\^  +  \E^\^) 


(62) 


The  expressions  for  Eg  and  E^j,  were  given  in  Eq.  47.  They  are  repeated  below: 


cot6> 

Eg  = - .  1  ^ - tan 


-jkor 


E<p  =  - 


^qDtMs  V  2s  y  2nr 

7’/o?7oMr(fy  cotOe~j'^°^ 


(63) 


D 


TE, 


2nr 


Combining  Eqs.  62  and  63,  we  have 
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The  last  integration  over  6  cannot  be  done  analytically.  The  directivity  is 


Furthermore,  we  can  define  the  radiation  resistance  in  the  usual  way: 


(67) 

We  would  also  like  to  calculate  the  power  launched  into  any  excited  surface  waves.  If  only  1 
surface  wave  is  excited,  then  it  is  feasible  to  use  the  expressions  for  the  surface  wave  fields 
derived  in  the  previous  section  to  determine  the  power.  However,  a  more  general  procedure  that 
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accounts  for  any  power  launched  into  surface  waves  is  used  here.  Consider  the  Poynting  theorem 
for  lossless  media: 


E  -r  dV^jj>  Savg  -dS+j^ III -  8\E\^)  dV 

If  we  take  the  real  part  of  both  sides  of  Eq.  68,  then  we  obtain 


(68) 


total 


1  Cff 

^  (fh  Savg  ■  ds  ^ --Re  E-J*dV 
JJ  ^  JJJ 


(69) 


In  Eq.  69,  we  take  S^vg  to  be  the  total  power  flow  produced  by  the  source.  Thus,  S^vg  in  Eq.  69 
is  not  the  same  S^vg  that  was  used  in  Eq.  64  to  calculate  the  power  radiated  into  the  upper 
hemisphere.  Working  with  Eq.  69, 

Ptotai  =  -^Re  lllE-(moS(x)S(y)S(z))dV 

00 

Ptotai  ^  jj  I  z=o  ^IoS(x)5(y)dxdy  (70) 

—  00 

By  Parseval’s  theorem,  Eq.  70  is  equivalent  to 

00  I 

rr  \ 

(71) 

In  order  to  use  Eq.  71,  we  have  to  detour  slightly  to  derive  a  convenient  expression  for  Ey..  Ey- 
may  be  found  from  the  vector  potential  as  follows: 


\  2=0 


dA 


jPx 


z=0‘’ 


(72) 


Erom  Eq.  29,  we  have 


—  n 


= - -e-ffcziz 


A 


Dte 


(73) 


Zi 


DxeD 


g  jkzyZ 


TM 


Combining  Eqs.  72  and  73,  we  have 

Ex  I  2=0  ~  ^^oVo 


Px  Pzy(^Er^r 

iEr{Pzy^  +  Py^y 

PqDteDtm 

PqDte 

(74) 


Equation  74  can  be  manipulated  into  a  form  more  convenient  for  the  integration  that  follows.  We 
begin  by  substituting  Eq.  34  into  Eq.  74,  yielding 
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^'^°[ik,^+k/)k,Dr,DrM  . 

From  the  definition  of  Dji^,  we  have 

^r^zi  =  -j'i^TM  +  K2  tan(/c22/i)) 

Combining  Eqs.  75  and  76, 

~|  .  ^Zi  [^Z2  "f  i/^r^Zi  (^TM  ”f  ^Z2  ^an(^/C22^))]  }kr(j^z^  +  ^y  )l 


z=0  =  ^^oVo 


{kj(  +  ky  ^koDf^Dj 


pf  ^Zi  [^Z2  J kr^Zi^i 

,2  tan{k^^h)  +  jiirkz^DjM]  jkr{kz^^  +  /Cy^)' 

.=0  -  ^m'/o  1^ 

~  1  l^x^^z 

p  —  />/  M  ^  ^iL 

k  ^  +  k 

2  ,  . 
•Z2  +]k 

^/Cgf^T’fi'f^T’M  ^O^TE 

r^Zi^Z2  ^an(^/C22^}]  ^  j kr^x  ^z^ 

\  z=o  ~ 

7Mr  (^Zi 

r  F  ky 

+fc/)l 

^k^DrE^TM 

(/C;^  +  ky  )/fo^T£ 

^0^T£  J 

<>f  ^Zi  [^Z2  "fj/^r^zi^ 

-Z2  tan(/Cz2/i) 

1]  cot(A:z2/i)  ^  ikrkx^K^^ 

-“-“"""1  (^  + 
jk-r(j^Zi  ”f  ^y 

ky^)koL 

!)1 

^TE^TM  COt(/ 

Cz2^)  +  ky^)koDTE 

k^Drs 

PI  ^Zi^Z2  [^Z2 

J 

^^/Cz2^)  j kr^z-il  1 

jkr^x  ^Zi  jkri^^Zi  ”f  ^y  ) 

°  °  L  +  ky^)koDTEDTM  CO 

^Zi^Z2[^T£] 

t(/Cz2/l)  (- 

1 

+  ky^)koDjE  koD-TE 

jflpkx  k^^  jk-rij^z-i  ”f  ^y  ) 

^■*0^0  /"  2  2  A 

Lv^x  +  ^y  J^oDte^tm  CO 
[/fx^fczifcz2  tan(/c22/i) 

t(/Cz2/l)  (/f 

,  ikrkx^K^^ 

+  ky^)koDTE  ^oDte 

—  jkri^X^  +  /Cy^)(A:z^^  +  /fy^) 

°^nfco(fc/+VKM 

_  „  |■/C;,^/Cz^/Cz2tan(/Cz2/^) 

r 

,  ikrkx^K^^ 

(/C;|.  +  ky  ^koDfE 

—  jkri^x^  +  ky^^ijc^:^  +  /Cy^) 

™[fco(fc/  +  VKM 

^>r  „  ^Zi^Z2  tan(^/C22^)  ^  jkrij^X  ^Zi  ~  ‘ 

(/c^f  +  ky  ^k^DjE 

^Zi  “  ^Zi  “  ^x  ^y  ~  ^y  ) 

^^°"n/co(fc/+VK 

~i 

Z7  —  pi  y.  A  . 

'M 

j^/Cz^  tan{k^^h) 

(^kx  +  ky  ^koDfE 

•ky^ipz^  +  ^x^  +  ^y^) 

~  1 

/2  1  —  y^/  Yi  _ 

'x^  +  ky‘ 

^k  k 

X  «-z^«-z2 

)^TM 

tan(A:z2h) 

(/c^f  +  ky  ^k^DjE 

jkrky  (^0  ) 

1  z=0  “  ^^oVo  ^ 

+  ^y  )^TM 

(ky^  +  ky  ^k^DjE. 
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(77) 


ko^kx^  +  ky^^Djj^  {kx^  +  ky^^DjEl 

Substitution  of  Eq.  77  into  Eq.  71  gives 

GO 

II 


total  ~ 


8n^ 


Re 


kx^k^^k^^  tan{k^^h) 

j ^0 

k^lkx  +  ky 

{k-x^  +  ky^^DjE 

d-  li  ^  d-  Ji  y 


(78) 


We  then  introduce  the  change  of  variables  into  Eq.  78: 

/Cv  =  kt  cos  a 


ky  —  kf-  sin  a 


(79) 


yielding 


total 


8n^ 


Re 


2n  00 

II 

0  0 


cos^  a  k^Jiz^  tSLulk^^h) 

Jlirko  sin^  a 

koDrM 

1 

kf-dkf-da 


total  ~ 


8n 


Re 


I 


kz-^kz2  jd-rko 


koD 


TM 


D 


TE 


kf-dkf- 


(80) 


In  the  limit  of  a  lossless  substrate,  poles  exist  along  the  integration  interval  of  Eq.  80,  thus  this 
integral  must  be  interpreted  in  terms  of  Cauchy  principal  value  integrals.  This  is  similar  to  what 
was  done  earlier  in  the  analysis  of  the  surface  wave  fields.  If  we  interpret  Eq.  80  as  a  contour 
integral  along  the  real  axis  of  the  complex  k^  plane,  and  we  adopt  the  same  definition  for  k^^, 

then  we  may  integrate  along  any  path  so  long  as  we  start  and  stop  at  the  same  points,  since  all  the 
functions  in  the  integrand  are  anal3^ic.  Consider  the  contour  shown  in  Eig.  4,  where  the 
integration  path  avoids  the  difficulties  associated  with  the  poles  by  integrating  around  them."^ 


Fig.  4  Alternative  path  of  integration 
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We  have 


Ptotai  -  Re(/i  +  l2+h+  U) 

However,  the  integrand  for  1 4  is  pure  imaginary  so  long  as  the  starting  point  is  past  so 

that  the  total  power  can  be  found  by  numerically  integrating  around  the  box  +  I2  +  I3.  The 
difference  between  the  total  power  and  the  radiated  power  is  the  power  launched  into  surface 
waves.  We  define 


Additionally,  it  is  worth  noting  that  integrating  from  0  to  /cq  gives  the  power  radiated  into  the 
upper  hemisphere.^  Although  a  proof  of  this  fact  is  not  known  to  the  author,  numerical 
calculations  were  always  in  agreement  with  results  obtained  from  Eq.  65.  Lastly,  we  may  define 
a  surface  wave  resistance  and  efficiency  as 

(82) 

(83) 

Some  numerical  results  obtained  from  the  formulas  in  this  section  are  included  in  Figs.  5-16. 
Note  that  wavelengths  on  the  horizontal  axis  are  measured  as  wavelengths  in  that  material,  and 
not  in  terms  of  free  space  wavelength.  In  all  cases,  a  dipole  whose  length  is  1/50  of  the  free 
space  wavelength  is  assumed. 

For  Sj.  =  1.01;  =  1.01;  £  =  do/50 
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Equivalent  resistances  of  horizontal  infinitesimal  dipole 


Height  h  (wavelengths) 


Fig.  5  Equivalent  resistances  of  a  horizontal  infinitesimal  dipole  situated 
on  a  dielectric  substrate;  =  1.01;  =  1.01;  ■£  —  Aq/SO 

Directivity  of  horizontal  infinitesimal  dipole 


0  0.04  0.08  0.12  0.16  0.2  0.24 

Height  h  (wavelengths) 

Fig.  6  Directivity  of  a  horizontal  dipole  situated  on  a  dielectric  substrate; 

£y  —  1.01;  Hr  —  1.01;  £  —  Aq/SO 
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Efficiency  of  horizontal  infinitesimal  dipole 
situated  on  dielectric  substrate 


Fig.  7  Efficiency  of  a  horizontal  infinitesimal  dipole  situated  on  a  dielectric  substrate; 
—  1.01;  =  1.01;  £  — 

For  =  10;  iiy  =  1.01;  ■£  —  Xq/SQ 
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Fig.  8  Equivalent  resistances  of  a  horizontal  infinitesimal  dipole  situated  on  a 
dielectric  substrate;  =  10;  =  1.01;  ■£  — 
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Directivity  of  horizontal  infinitesimal  dipole 
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Fig.  9  Directivity  of  a  horizontal  infinitesimal  dipole  situated  on  a  dielectric 
substrate;  =  10;  fij.  —  1.01;  ^ 
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Fig.  10  Efficiency  of  a  horizontal  infinitesimal  dipole  situated  on  a  dielectric 
substrate;  Ej.  —  10;  fij.  —  1.01; 

For  Sj.  =  1.01;  [ij.  =  10;  £  =  Xq/SQ 
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Equivalent  resistances  of  horizontal  infinitesimal  dipole 


Height  h  (wavelengths) 


Fig.  1 1  Equivalent  resistances  of  a  horizontal  infinitesimal  dipole  situated  on 
a  dielectric  substrate;  =  1.01;  =  10;  ■£  —  Aq/SO 
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Fig.  12  Directivity  of  a  horizontal  infinitesimal  dipole  situated  on  a  dielectric 
substrate;  —  1.01;  —  10;  f  —  Aq/SO 
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Efficiency  of  horizontal  infinitesimal  dipole 


Height  h  (wavelengths) 


Fig.  13  Efficiency  of  a  horizontal  infinitesimal  dipole  situated  on  a  dielectric 
substrate;  =  1.01;  =  10;  f  —  Aq/SO 


For  =  10;  /Tj.  =  10;  ■£  =  Xq/SQ 

Equivalent  resistances  of  horizontal  infinitesimal  dipole 
situated  on  dielectric  substrate 
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Fig.  14  Equivalent  resistances  of  a  horizontal  infinitesimal  dipole  situated  on 
a  dielectric  substrate;  =  10; =  10;^  =  Ao/50 
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Directivity  of  horizontal  infinitesimal  dipole 
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Fig.  15  Directivity  of  a  horizontal  infinitesimal  dipole  situated  on  a  dielectric 
substrate;  =  10;  fij.  —  10;  £  —  Aq/SO 

Efficiency  of  horizontal  infinitesimal  dipole 
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Fig.  16  Efficiency  of  a  horizontal  infinitesimal  dipole  situated  on  a  dielectric 
substrate;  Ej.  —  10;  iij.  —  10;  =  /^o/SO 
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7.  Conclusion 


In  this  report,  the  solution  of  the  boundary  value  problem  associated  with  dipole  radiation  over 
grounded  magnetodielectric  substrates  was  treated.  The  mathematical  details,  which  proved 
difficult  to  find  in  the  literature,  were  reproduced.  The  analysis  of  these  solutions  in  the 
subsequent  sections  derived  the  radiation  field,  surface  wave  fields,  directivity,  radiation 
resistance,  and  efficiency.  Some  numerical  results  were  plotted  in  Figs.  5-16.  From  this  small 
sampling  of  data,  we  can  already  conclude  that  the  behavior  of  planar  antennas  printed  on 
magnetodielectric  substrates  is  very  complicated  and  non-intuitive.  For  example.  Fig.  14 
indicates  that  a  very  small  dipole  printed  on  a  very  thin  substrate  could  easily  present  an  input 
impedance  of  50  ohms  at  its  terminals.  However,  despite  conventional  wisdom  that  quarter- 
wavelength  substrates  make  very  effective  radiators.  Fig.  14  indicates  that  essentially  all  of  the 
available  power  is  delivered  to  surface  waves  excited  within  the  substrate  itself,  which  will  not 
necessarily  lead  to  desirable  performance. 

Furthermore,  a  novel  treatment  of  the  complex  valued  square  roots  was  shown  to  produce  results 
in  agreement  with  those  discussed  in  the  literature  with  an  arguably  less  tortured  mathematical 
procedure.  Appendix  B  can  be  consulted  for  an  explanation  of  this  approach  and  a  derivation  of 
the  properties  of  the  complex  valued  square  roots  in  the  complex  plane.  This  should  find 
application  to  other  boundary  value  problems  (such  as  the  present  one)  for  which  solutions  are 
obtained  using  Fourier  transforms  and  those  that  require  the  determination  of  inverse  transforms 
to  extract  useful  formulas. 
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The  boundary  conditions  used  in  the  derivation  of  the  magnetic  vector  potential  are  developed  in 
this  appendix.  The  preliminary  assumption  made  in  the  analysis  of  the  problem  was  that  the 
electromagnetic  field  could  be  modeled  by  2  components  of  the  vector  magnetic  potential: 


A  —  xAx  +  zA^ 

The  electromagnetic  field  associated  with  this  vector  potential  is 

1 

H 

1 

E  = - V  X  H 

j(x)£ 


Substituting  Eq.  A-1  into  Eq.  A-2  yields 


1 

H  =  - 


X- 


+  y 


/dAy- 

V  dz  dx  ) 


dA, 


E  = 


JO)  [IS 


jLl  dy  V  dz  dx  J  dy  J 

_  d^A^  d^A;^  d^Ax\  ^/d^A^  d^Ax\ 

l^dxdz  dz^  dy^  J  ^  ygyg^  gxdy J 


+  z 


f 


d'^Ar  d^A,  d'^A, 


\dxdz  dx^  dy^ 

Given  the  general  form  of  solution,  Eq.  A-3  may  be  reformulated  as  follows: 


1 

H  =  - 


-x{]kyA^)  +  y  (^+jkxA^^  +  z{jkyA^) 


jOOJLE 


dA 


dA, 


X  (  -jk^  +  (/c^  -  k^^)A^  )  -  y  ( jky  +  kJiyA 


dz 

\  / 

+  z  (-jkx  ^  +  {k^  -  k^^)A'^ 


(A-1) 


(A-2) 


(A-3) 


(A-4) 


(A-5) 


At  =  — /i ,  a  PEC  interface  requires 


zxE  =  0 
z-H  =  0 


Substituting  Eq.  A-4  into  Eq.  A-7,  we  immediately  obtain 


z=-h 


=  0 


Then,  if  we  apply  Eq.  A-6  to  the  y  component  of  the  electric  field  in  Eq.  A-5 


dA.  _ ' 

_ fi  —  k  k  A 

g^ 


=  0 


(A-6) 

(A-7) 


(A-8) 


(A-9) 
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Use  of  Eq.  A-8  in  Eq.  A-9  yields  the  other  boundary  condition  applicable  at  the  PEC  interface: 


dA 


Z2 


dz 


z=-h  ~  0 


Now,  we  consider  the  boundary  condition: 

=Js 

In  reality,  there  is  no  surface  current,  yet  reconsider  the  current  density: 

/  =  x£lQ5(r)  —  x-£lQ5(x)S(y)5(z). 
Using  the  integral  representation  of  the  delta  functions 


—  00 


e-j’^y^dk. 


we  have 


(2tt)‘ 


We  identify 


and  therefore 


uu 

jj  [2-f/od(z)]e“-'^  e~^^y^dky.dky 


J  =  x£Io5(z)  =Js5(z) 


z  X  (i/i  -  H2)  =  x-eio 
Using  Eq.  A-4  in  Eq.  A- 16  yields  2  equations: 


A 


Zi  I  z=0‘' 


=  —A 
kr 


Z2  \  z=0 


'dA. 


1  fdA 


jky'A^  z  =  0  ..  j^X^Zo  z=0  ~~ 


y  dz 

Rearranging  Eq.  A-18, 

/  dAx,  1  dA 


fij-  \  dz 


X2 


y  dz  fij-  dz 
Using  Eq.  A- 17  in  Eq.  A- 19,  we  obtain 


^  I  Z  =  0  j^X  (^^Zi  ^  ■'^22^  I  Z=0  ~  ko'^^0 


(A-10) 

(A-11) 

(A-12) 

(A-13) 

(A-14) 

(A-15) 

(A-16) 

(A-17) 

(A-18) 

(A-19) 
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dA, 


1  dA 


X2 


dz 


z=0'’ 


fly  dz 


2=0-  “ 


(A-20) 


We  can  obtain  another  boundary  condition  on  the  vector  potentials  from 


f^r 


z=0 


=  0 


which  gives 


=  A, 


z=0 


Our  final  boundary  condition  must  come  from 

2x(£i-®2)I„„  =  0 

Applying  Eq.  A- 23  to  the  y  component  of  the  electric  field  in  Eq.  A-5,  we  find 


(A-21) 


(A-22) 


(A-23) 


-jK 


dA 


Zi 


Y  ^  dz 
Rearranging  Eq.  A-24, 


kxkyAx^  1  I  2=0 


jlyEy 


-jky 


dA 


Z2 


dz 


kxkyAx2  1  I  2=0 


=  0 


(A-24) 


-jky 


'dA 


Zi 


dz  pLfEy 


A. 


dA 


Z2 


dz 


z=o 


Substitution  of  Eq.  A-22  into  Eq.  A-25  gives 


dA 


Zi 


1  dA 


dz  /Ll^Er  9z 


Z2  I  _  kr^r  1  j 

,_Q-  —  Jf^y  -^1 


Ey 


1  I  z=0 


(A-25) 


(A-26) 


In  total,  there  are  6  unique  boundary  conditions  applicable  to  the  magnetic  vector  potential.  They 
have  been  derived  in  this  appendix  from  the  boundary  conditions  satisfied  by  the  electric  and 
magnetic  fields.  They  are  repeated  together  as  Eq.  A- 27. 


(A-27) 
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The  definition  of  for  real  values  of  kf-  is 


,  kf-  <  kf) 


(B-1) 


In  order  to  employ  the  methods  of  complex  integration,  this  definition  must  be  extended  to  the 
whole  complex  kf-  plane.  The  definition  we  use  for  k^^  is 


kt^\exp  |argo(/co' 


(B-2) 


In  Eq.  B-2,  argo(z)  is  the  value  of  arg(z)  that  lies  in  the  interval  (0,27r].  This  is  depicted 
graphically  as  follows: 


tt/2 


3tt/4 


TT  - 


n'/4 


5tt/4 


7tt/4 


3tt/2 


From  the  graph,  we  see  that  argo  (z)  is  discontinuous  if  we  cross  the  positive  real  axis  in  the  z 
plane.  The  real  axis  is  said  to  form  the  branch  cut  for  argo(z)  and  argo(z)  is  said  to  be  a 
particular  branch  of  the  multi-valued  arg(z)  function.  Now,  z  =  k^^  —  k^^,  so  that  we  may  map 
the  branch  cuts  from  the  z  plane  to  the  kf-  plane.  The  real  and  imaginary  parts  of  z,  k^,  kf-  are 
denoted  as 


z  —  X  +  jy 

^0  =  ■*"  j^Of  (B-3) 
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Then, 

x+j-y  ^  +;2(V^0i  -  (B-4) 

The  positive  real  axis  in  the  z  plane  corresponds  to  the  conditions 


z  =  X  +  jy 
X  >  0 
y  =  0 

such  that  the  equations  defining  the  branch  cuts  in  the  kf-  plane  are  seen  to  be 


(B-5) 


(B-6) 


The  second  of  these  equations  defines  a  hyperbola  in  the  kf-  plane: 


•'I 


(B-7) 


Since  kQ.  must  be  a  negative  number  to  give  attenuation,  this  hyperbola  will  appear  as  shown: 


From  Eq.  B-6,  the  real  part  must  be  greater  than  0.  The  endpoint  of  our  hyperbolic  curves  must 
occur  where  z  =  0.  This  corresponds  to  the  points  k^  =  ±/co-  Our  branch  cuts  are  as  follows: 
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In  the  limit  of  no  loss,  kQ.  0.  Equation  B-7  indicates  that  in  this  limit,  the  hyperbola  will 
shrink  inward  to  the  coordinate  axes.  The  governing  equations  for  the  branch  cut  reduce  to 


^Or  ^tr  +  —  0 

=  0 


(B-8) 


These  equations  admit  2  solutions.  That  is. 


kf-^  =  0  (the  imaginary  axis) 

kf-.  =  0  ]  \kt^\  <  \ko\  (the  line  segment  on  the  real  axis  from  —  k^  to  +  k^) 

This  is  consistent  with  the  physical  picture  of  hyperbolas  collapsing  inward  toward  the 
coordinate  axes.  The  branch  cuts  are  depicted  as  follows: 


(B-9) 


Branch  Cut  — 

In  view  of  Eq.  B-1,  and  in  order  to  obtain  convergent  integrals,  we  insist  that  Im(/c2^)  <  0.  We 
must  verify  that  definition  Eq.  B-2  satisfies  this  requirement.  In  quadrant  1,  we  know  that  k^-^  > 
0  and  kf-.  >0.  So  that. 


z  =  x+jy^  -j'2-kt^kt^ 

y  <0 


Thus,  in  quadrant  1,  we  have 


and 


I  <  ^argo(/co^  -  ^  tt 


jexp 

|argo(/co^  -kt^)  j 

[exp 

^argo(/co^  -kt^)  ] 

(B-10) 


(B-11) 


(B-12) 
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In  view  of  Eq.  B-2  and  Eq.  B-12,  we  can  conclude  that  within  quadrant  1,  the  following 
inequalities  hold: 

Re(/CzJ  >  0 

I 

<  0 

This  analysis  can  be  repeated  for  quadrants  2,  3, 4.  The  topology  for  is  summarized  as 
follows: 


Re(fc,J  <  0 
lm(/c^J  <  0 


Re{k,J  >  0 
<  0 


k 


0 


Re(/c,J  >  0 
lm{k^J  <  0 


Re(/c,J  <  0 
lm(fc^J  <  0 


Branch  Cut  — 


-13) 
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